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Preface 


This report documents the progress on the research being performed 
by the School of Aeronautics and Astronautics, Purdue University, for the 
NASA Ames Research Center/Qryden Flight Research Facility, under Cooperative 
Agreement/Grant Number NCC2-288*. The objective of the program is to 
determine the appl icabil ity of several multi -input/multi -output control 
synthesis techniques for the synthesis of flight control lav/s for advanced 
experimental aircraft, and to extend the techniques as necessary to obtain 
simple, robust control laws that meet specific handling qualities objectives. 

This report includes results from the synthesis of control laws for 
an advanced STOL vehicle in a low-speed approach flight condition. Two of 
four candidate synthesis techniques are reported herein - direct eigenspace 
assignment and explicit model following via a Linear Quadratic Regulator 
(LQR) formulation. The other two techniques under current investigation 
are implicit model following (also via LQR) and Cooperative Control 
Synthesis. This latter technique has been a topic of continued research at 
Purdue University, and utilizes pilot-in-the-loop (pilot modeling) 
techniques. 

A fundamental objective in this work is to obtain low-order feedback 
compensators, synthesized- via the techniques above, and with judicious use 
of state-estimation thus allowing the use of a reasonable number of sensors 
for feedback. This report focuses on the results from the synthesis of 
state-feedback control laws. However, two state-estimation techniques 
under consideration are noted. In implementing these estimation techniques, 
we are interested in the ultimate stability robustness of the system, and not 

increasing the dynamic order of the stick-response transfer functions: 

★ 

Mr. E.L. Duke, Technical Monitor 



The latter goal is motivated by the desire to preserve the handling 
characteristics "built into" the state-feedback control law utilized. 
Results on the implementations with state estimation will be documented 
in the next progress report. 



Contents 


page 


1. Introduction 1 

2. The Study Vehicle 3 

3. Direct Eigenspace Assignment 8 

4. Explicit Model Following 30 

5. Implementation Using Output- Feedback 46 

6. Conclusion 47 

7. References 49 

8. Appendix - MATRIX^ Macros 51 



1 


i- INTRODUCTION 

The ultimate performance objective of the flight control design is 
to make the combined pilot and vehicle system behave suitably. This 
objective introduces the idea of "handling quantities"* defined as 
those desirab le_ dynami c traits of an aircraft that wi.ll help the pilot 
and vehicle system perform their intended mission. With the "handling 
qualities" specifications available to the control system designer* he 
must use these specif ications to formulate and implement the control 
system design* and finally to verify - that- the- augmented- vehic le does 
- indeed meet: the-hand 1 ing .. qual i t i es. specifications. 

In the past decade or so* many neu> methodologies for - designing con- 
trol systems for multi-input multi-output systems have ;emerged. The 
objective of this study is to investigate the applicability of some of 
these 'modern' techniques to design of flight control systems* with the 
specific objective of meeting the handling qualities reqirements. 

In this particular study* tuo of these techniques — direct eigen- 
space assignment (DEA) and explicit model following (EMF)* are used ini- 
tially to synthesize control laws for the logitudinal dynamics model of 
a STOL vehicle in the landing approach configuration. The vehicle model 
and the flight control design requirements are presented in Section 2. 

In the succeeding sections* the two synthesis techniques are briefly 
discussed and the handling qualities specifications mapped into the 
~ algorithm formulation — The -control laws resulting from exercising the 
algorithms are evaluated in terms of achieved performance and robust- 
ness. 
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Since the synthesized control laws involve full-state feedback* 
methodologies for implementing' these' control 'lauis-using output-feedback 
and uithout adversely affecting performance and robustness* are a topic 
of significant interest and are presently being pursued. Promising 
techniques being considered are briefly discussed in Section 5. 

Finally* the salient features of the two design techniques are sum 
marized and the areas that require further investigation are suggested. 
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2. THE STUDY VEHICLE 

2 . 1. Vehicle Model Description ' 

The vehicle considered is a STOL aircraft, with an airframe similar 
to an F— 18 aircraft. The linearized dynamic model for the longitudinal 
axis includes four rigid-body degrees of freedom and three first order 
actuator lags (each at 15 rads/sec). The trim values for the landing 
approach flight coridition are listed in- Table- 2. 1. The control inputs 
to be utilized are the horizontal tail (elevator), the thrust vector 
angle for a 2-D nozzle, and the trailing edge flap. The state vector for 
the vehicle model is 

*r% Sh; Stv> Sf] 


where 


u.1 

U-fi = nondimensional forward speed ( -jj- ) 

K = angle of attack (rads) 

— pitch rate (rads/sec) 

9 =* pitch -angle (rads) 

= horizontal tail deflection (rads) 

&tv = Thrust vector angle (rads) 

Sf = trailing edge flap deflection (rads) 


with the control inputs taken as commands to the servo actuators, or 

S c ' [&Hc , ^' tV c , ^Fe 1 

The vehicle responses of interest" are 


where 
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fig. = normal acceleration at C. C. <g 's# +ve up) 
T = flight path angle (rads) 


The vehicle dynamics can then be written in the form 

X - fix +& Sc. 
f-C* 

The. open loop. system matrices are listed in Table 2.2. The open-loop 
vehicle eigenvalues are 

0 \557 j: j 0*1962 , -1-53 5 , = 0-8532 

Note the unstable pole. 


2. *2 Flight Control Design Requirements 

The flight control design requirements can be summarized as^ follows 

* Explicitly include handling quality criteria 

* Avoid excessive control surface rates 

* Reduce control energy at high frequencies 

The handling quality criteria is stated in terms of desired short 
period pitch rate response and the phase relationship between the flight 


path angle and the pitch angle i. e 

<v (%\ 

F' \ Fa 






) 


and 


r 

9 



with 


{% ^ - o-5 

V Q JjLh * 


= o-7s .mu 


J 
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and h). L and 

requirement C1D and a rise time for pitch rate# 0.8 sec. Furthei — 

more# the augmented vehicle should exhibit classical, “phugoid" like 

i 

dynamics for the lout frequency mode. 

In the present study# our primary objective will be to meet the 
above handling qualities specifications through proper control augmenta 
tion of the vehicle. 


selected so as to satisfy Level 1- handling qualities 




Table 2 . 2 Open Loop System Matrices 


7 


CM »-« «-• O 

MhM O 

noo o 

O ^ CM o 

o d o o' 6 6 ui 

ii ^ 


in n in o 

CD co cm o 

o«tn o 
o o «a- o 

o’ o’ h o o’ in o' 

ill 

i 


*-• O N 

in *1 cm 

CO O' *0 

O O N 


o o r-« o in o o 

I I 1H 


0 * O r* 
00 PI >0 

in h o 

y* O O 


0 o o o o o o 

1 I 


H (U rt o 

CM in CM O 
pi »-« o 
«MO 


000^000 
I I 


P) N PJ 
CM PJ O' 
O CM O 
CM PJ ^ 


00^0000 
I I 


in o oo 

NOO 

in o o 
o cm >o 


0 o o o o o o 

1 I I 




0000 



8 


2- DIRECT EIGENSPACE ASSIGNMENT - (PEA) 

3. i Gain Sunthes is 

In the direct eigenspace assignment technique* the control objec- 
tives are stated in terms of a desired eigenstruc ture for the augmented 
system. For the state feedback case* the synthesis problem is as follows 


Given a system 




system dynamics) 


u. = u c +Up j . a e fi 


7a 


(total control input) 


Up <= pilot's control inputs 

cf " 7 fy' ^ ^ (system responses) 


find a control law of the form 




to achieve some desired eigenspace for the augmented system 

X = (A - 6 1 i) x + 6 Up 

y-a* 

To determine K* note that the augmented ( c losed— loop ) system eigen- 
values and eigenvectors are related by 

(A-Bk) y H -- > c . L--I,., n 


(31 -l) 


where 

^ = closed-loop eigenvalue 
>£.°* closed-loop eigenvector 

For full-state feedback* the limitation on the achievable eigen- 
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space is that all the 7 desired' closed-loop - eigenvalues can be exactly 
placed while only "m" elements of their associated eigenvectors may be 
specified (where m = dimension of ED. Since in general m < n» we cannot 
exactly obtain all elements of the desired eigenvector for each closed- 
loop mode. One approach is to determine the " bes t " . ach i evab le eigenvec- 
tor * for each of the closed— loop modes* that minimizes the mode's 
1 , 

cost function 




C -)}•< • y ^ 




where 


= i'th achievable eigenvector associated with eigenvalue 
^ = i/th . desired eigenvector 

= i/th n-by-n- symmetric -.posit ive- semi-definite weigh ting 
matrix on eigenvector error elements 


and * denotes conjugate transpose. 

Equation (3. 1. 1) can be rewritten as 

Defining the vector Uf- ~ -fcy and using eqn. (3.1.3)* the solution to 

. a i : 

(3.1.2) is obtained as 

U ti-Tftrfi t-i] 

where £>. 


□nee are obtained* the achievable eigenvectors are given by 

I* 

\= Li'*i . i’<>- >* 

and the feedback gains are obtained as 

K'-WV ' 1 
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where 

W = [ &, ^ J 

V = [fy, •; • • ,^y» 1 

This solution algorithm is due to Schmidt and Davidson C21. A 
macro using MATRIX.. commands uias written to implement the above syn— 

t 

thesis procedure. This macro is documented in the Appendix. The block 
diagram for DEA control law implementation is shown in Fig. 3.1. 


3 • 2 Selection of Desired Eiaenspace and Resulting Control Law 

The control design objective is to make the augmented vehicle modes 
'like those of a- “classical "_airp lane. i. e... the phugoid and short period 
modes.' For this- initial- synthesis* we~ chose- the-*desired phugoid mode to 


D 0 — 

= o- Z5 McU /a&c , 

^ = -o-odtj. 0-^87 

with the corresponding desired eigenvector selected as 

VdJ - o , a , a.±^l , a. , a] [turi* Jvf-, p]- £) 

where "a" denotes an arbitrary value. The' above choice for the phugoid 
eigenvector reflects the desire that the phugoid mode shape be dominated 
by forward speed and pitch attitude response with little or no angle of 
attack contribution C41. 

The short-period mode frequency and-dampingare selectedto reflect 
the handling qualities requirements. This choice is as follows 

LOa^ - 3-7 Mu h> /Me t = o • 7S~ 

=> ^ =' ' - a-7% i $.3^3 

The choice of the desired short-period eigenvector is based on the 
requirement that the short-period mode be dominant in angle of attack 
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and pitch rate response* with little or no forward speed contribution. 
Furthermore* the short-period eigenvector should reflect the desired 
flight path angle to pitch attitude phase relationship i. e. we want 



Using the relationship f=9'0(> we can write the above constraint as 



With O'lS' Mk. and ~7‘7& - ^ S-S'i , we get the desired relationship 

between 

0( and in the short-period mode to be ( 

= - 0-181 + j, 0-3 o38 

Then an eigenvector -that results is 

hT •' [ 0) ~0l%\+ j. 0'3o3%j I,. a J '(toil sb-k pfr?) 

The above choice as well as the desired eigenspace for the actuator 
modes is listed in Table 3.1. The actuator poles are left near their 
open loop values* and the corresponding eigenvectors are selected simply 
to obtain decoupled actuator modes. 

From Table 3. 1 we note that for the short period mode and for all 
the three, actuator modes, we are in effect specifying only three ele— 

- ments of the corresponding eigenvectors. Since we have three controls 
.. available* we shall be able to exactly achieve the specified elements of 
these eigenvectors. The phugoid eigenvector effectively has only two 
specified elements. This leaves one extra degree of freedom to further 
constrain the. phugoid eigenvector. One possibility to be explored in 
the future is to use this freedom to reduce the control surface deflec- 


_ • >v 

* . <*> ,.'4 * 



tions. 





The achieved eigenvectors corresponding to the above selections are 
listed in Tab le. 3. 2*__and the assoc iated . feedback control gains in Table 
3.3. As we expected* the desired eigenspace is achieved exactly. For 
comparison* the phasor diagrams for the open— loop and achievable eigen- 
vectors are shown in Fig. 3. 2. These indicate "classical" phugoid and 

short-period mode like behaviour for the augmented airplane. Note* how- 

» 

ever, ^the large control deflections in these mode shapes* indicating 
high deflection requirements for this control law. 


3. 4 Performance Evaluation 


*» when calculated from these transfer functions for J, X does 


With the pilot's input taken to be commanded horizontal tail 
deflection' ( th e~ transfer functions between the responses of 

interest -and the.-pilot's input are listed in Table 3.4. These transfer 
functions indicate that the augmented system response exhibits the 
desired decoupling between the phugoid and the short period mode. More- 
over! 

give the value (-0. 181 + jO. 3038) that was specified during the design 

process. Also in the transfer function* note that has increasd 

f>c : -i 9 3 l 

to 0.92 from its open loop value of 0.5 sec . 

Other responses of interest are the normal acceleration at the 
center of rotation <^2 CR ) and the flight path angle at the center of 
rotation (~^c.r.). With the horizontal tail deflection as the only pilot 
input, the center of rotation is located 10.45 ft ahead of the C. O. The 
transfer functions for these responses are; also listed in Table 3. 4 

The transfer functions from pilot's control input to the actual 
control deflections are listed in Table 3.5. 

Time histories for a step pilot input (£^=1°) are shown in Figs. 
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3.3 to 3. 5. Fig. 3.3 shous 
- these again indicate the 
behaviour. Fig. 3. 4 shous 
flight path angle* both at 
Fig. 3. 5 shous the actual c 


the response for the augmented vehicle states 
decoupled phugoid and short-period mode 
the response for normal acceleration and 
the C. 0. as well as the center of rotation, 
ontrol deflections for this step pilot input. 


3.- & Robustness Evaluat ion 

A very important consideration in flight control design is the sta- 
bility robustness of the augmented system. Given a system uith transfer 
function matrix G(s>* uhere ^ 

G,CO B> 

a reliable (but sometimes conservative C53) measure of 
minimum singular value of the return difference matrix* 
function of jliJ * or in this case* uith the loop broken 

2T [ 4 + < <1 CO ] 

For the control lau obtained above* the singular values of the return 
difference matrix are plotted in Fig. 3.6. Note the relatively lou sta- 
bility margins in the frequency range of upto 4 rads/sec (the frequency 
range of the vehicle* except for the actuators). Hou to increase this 
robustness requires additional consideration.. 

A more "classical" approach to evaluating the stability robustness 
is to "break" one loop at a time (i.e. one loop open and the other loops 
closed). The eigenvalues of the system uith the various loops open* and 
for the above control lau* are listed in Table 3.6. (Note that the sys- 
tem uith loop open has one unstable pole* uhich means that the con— 

i . 

trol lau uill not be able to stabilize the vehicle in case of an eleva- 
tor actuator failure). The Bode plots for the one-loop-at-a-time 


robustness is the 
evaluated as a 
at the input* 


f 
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analysis are shown in Figs. 3. 7(a)-<c >. - These plots appear to indicate 
"good" gain and phase margins in all the three control loops. This 
might be optimistic in light of the stability margins based on the 
singular value analysis. Whether this controller has sufficient robust 
ness is an open question. If not< how to improve robustness is a topic 
of interest. 



>-w># 


15 






«v 




16 


u 

o 

4 * 
u_nj 
I/O 3 
4* 
U 
< 


(VI 


in 


ca <*■ o 
CM N *<■ O 
O O O 

o o o o 


I 

w 


o o o 
. I 


O O O th 


t- 

o 

cl 

u 

< 


if) 


n n n >o 
n o «o 
o o o o 
o o o o 


o o o o o -« 

i i 


u- 

oo 


ail 


v- 

vo 


r 

<-0 


T 3 

CL 

in 

c 

at 

on 


. x 
U 3 


f- 

O 

-P 

m 

3 

+> 

U 

< 


^ h n n 
nnino 
oonio 
0 0^0 


o o o 

i_LJ 


0 ^ 

1 


o o 


cb 


iu 


ai 

> 

ai 




01 

O 

ai 


•> 

cr" 

■a 


1 03 

n : 
' o 

cm ^i n 1 
oo CD O' n 

NfHif 


o 

•r < 

Cl 

*r 

n 

^ o N O 

^ . 

f- 


6 

deinn 


a 

oi 


n n -7 -> 

|3 

a. 


if 

h- H- -H H- 

l 1 


fi 

0* 

O' CD CD in 


+* 

03 

o 

n cn n m 

v* 

(- 

fs 

CD 

o in ^ in 


o 


H 

oi -< oi «r 

£ 

01 

• • 

• » • • • 

cn 

1 

0-0 ■*-* OOhOI 



w 

1 ! 

11 » l 


~ r 


(Q 

»- 


N 

as 

01 


Dsns 
N O CD O 
0* O' 0* N 
^CONN 


N 01 O' 

o «■ oi 
-i n o 

H Q* H 


■w o 


c o *+ o o <m m 


o 

o» 

3 


+1 

in 

o 


If If If -H 
<«hO 

0 o oi m 
in h >o s 

01 o ^ o 


H- +1 1 + 
*0 CO **■ 
MtO> 

0 0 * CJ 

01 N N 


W* O O N 
I I 


0 o o 

1 I I 


,?'/A 

' is ’ 



17 




1 N n o * 
n h cy 
<r «r *o 
o o o 



6 o' 6 

i 



(> cd n 

■ c 


o o cm 
000 

H O H 

- odd 

CJ 


1 

•H 



01 

0i 

“O 


N CM 00 

< 


O -rt 

in h 0 

- UJ 


•-=.- cm 0 cm 

Q 


• • • • 

t- 

o 


000 
1 1 

<4- 

u. 


01 

\^0 

H H 

c 


O N CD 

■H 

. H 

0 >0 >0 

<n 

V/O 

*-« 0 n 

CTl 

<■* 

• 


ar 

000 


oo 

1 , 1 

u 



hj 

XI 

<3} 


"O 


CM CM O' 

01 

o^- 

«o «r cd 

01 


>o ^ cn 

IL 


cm in in 



CM O 'H - 

CO 

id 

! 

co 

U-..-I 


01 

V-" 

r. O' cm 

H 


'O cm n* 

XI 


ti CD -0 

<o 


h n 0 

H 


• • • 



O in cm 



U '-t 1 



1 

CM CM -0 



a) CM CD 
CM CD CD 



^ CM CD 



H O H 



t 



Table 3. 4 Transfer Functions for DEA desi 
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Table 3. 5 Control deflection Transfer functions for PEA desi 
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Figure 3. 2 Eigenvector Phasor Diagrams 
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Figure 
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Figure 3.7(c) Bode plots for op loop open 
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4- explicit MODEL FOLLOWING ( EMF ) 

4- i Gain Sunthesis 

In the LG explicit model-fol lowing algorithm* the problem formula- 
tion is as follows 

Given vehicle dynamics 

tc v *-.A v %v +& v “ 

2 V - 

where 2"^ is the vector of vehicle responses to' be controlled* and 
desired system characteristics in the form of a state-space model 

" Ayn ’Z'y r \ •+• £*£ 

r C-m 


where <$^.= pilot stick input modelled as 

«.* * ' % ^ + % ** 

we want to determine UL to minimize 




t e = E fa if J T [^v - S* ^ + ^ “ ] dt ! 

T-9 00 o 

Combining the vehicle dynamics and the model dynamics* we get 


Y>v 


°At 


Av 0 

Am 
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X ■> 

y~7r\ 
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B v 

O 





_T r- "T "| 

Writing the augmented state vector as % ' L x v ) ^hw^atJand with appropriate 
definitions of the matrices A* B* C and D* we can write the combined 
system dynamics as 


x - f\x -v- 8 ul -l J) 

?- ■= Z ' ? 9 Cx 

Z-n t- .. £■ Tin 
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and 3^ can be written as 

T , E 4 f [x r c r Q^cx +r a T R«] 1 

e ‘ T->00 o' 

The problem is now in the form of a Linear Quadratic Regulator (LQR) and 
the resulting control law is C63 



with 

k=i cVp - 

and P^o and symmetric* the solution to the algebraic Ricatti equation 

/f P + PA + C T 6( ? C -Pe/T 1 6 t P = 0 

The block diagram for EMF control law implementation is shown in Fig. 

4. 1. A-macro- using- MATRIX j^-commands^for-the -synthesis of the EMF con- 
trol law is documented in the Appendix. 

The augmented system results in the following state-space represen- 



The vehicle state-space representation is as in Section 2. 1 with 
all the three inputs ( 6 h c an< * ve h* c l e responses to be con- 
trolled. J . 


With the handling quality criteria stated in terms of desired short 
period response for the augmented system, we have 


■ V = y? ex )_ . 



3 Z 


u+ y? e± .'> 


In the state-space form* this model can be represented as the third 
order system 
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(based on stick gradient of t°/lb);- the numerical values 

of the model . system matrices are listed in Table 4. 1 

For the purposes of the following control law synthesis* the stick 
time constant 

and the output error weighting matrix Q = I were chosen. The 


control weighting 

matrix R 

was chosen 

as 
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where subscript "max" refers to the maximum allowable deflection. As 
the "control authority" is increased (scalar weighting $ is decreased)* 
the crossover frequency of the loop transfer functions k v ^-/^3 will 
increase. Higher loop cross-over* in general* relates to a higher 
fidelity match of the vehicle and model. However* higher crossover fre- 
quency can also increase the chance of the control system exciting the 
unmodelled modes (i. e. structural modes) and. may lead to undesirably 
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high control deflections. A value of S = 1 was chosen for this prelim- 
inary investigation. Although the resulting loop transfer functions are 
not presented here* the cross-over frequencies for the chosen value of f 
were found to be quite (* fO rad/s) 

For the above problem formulation* the feedback gains obtained by 

exercising the solution algorithm are listed in Table 4.2. 

\ 

4 . 3 Performance Evaluat i on 

For the control law synthesixed as above* the transfer functions of 

V 

interest for the augmented system are listed in Table 4.3. From the -^r 
and -r- transfer functions* we note that unlike the DEA case. model reduc- 
tion is required in order to get the low-order equivalent system parame- 
ters to determine the handling qualities of the augmented system. 

The desired frequency responses (from the model) and those for the 
augmented vehicle are compared in Figs. 4.2 to 4.4. Fig. 4.2 shows that 
the pitch rate frequency response for the vehicle closely approximates 
the model response for frequencies upto 10 rads/sec. Fig. 4.3 shows a 
fair agreement between the desired flight path angle response and that 
obtained at the center of rotation through augmentation. Fig. 4.4 com- 
pares the resulting flight path angle to pitch attitude relationship 
with that desired. This last result is quite unsatisfactory and a 
redesign with higher control authority may be required in order to 
obtain a better fit between the model and the augmented vehicle. 

Time histories for a step pilot input (<5^. = /°) are shown in Figs. 

4.5 to 4.7. From the 0( } ^ and S responses in Fig. 4.5* we note that 
the augmented vehicle does exhibit "classical" aircraft like dynamics. 

The normal acceleration and the flight path angle response* both at the 
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C. G. as well as the C. R. # are shown in Pig. 4.6. The control deflec- 
tions for a step stick input are shown in Fig. 4.7. 

4- 4 Robustness Evaluation 

Using the block diagram of Fig. 4. 1» the stability robustness of 
the EMF design may be evaluated by considering the values of £3+ kv ^ (j-^l 
where G) v Ci) = [*6 6 V . The singular values of this return difference 
matrix for the design obtained above are shown in Fig'. 4.8. Note that 
these stability margins are better than those obtained for the DEA 
design. If even higher stability margins are desired* then a control 
redesign with an appropriate choice of the control weighting matrix R 
may be required. 



Table 4. 1 Model System Matrices 
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Table 4. 5 Feedback gains for EMF design 
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Table 4. 3 Transfer Functions for EHF desi 
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Figure 4.7 Control Deflections for step pilot input (fu'l*) 
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Figure 4.8 Singular values of return difference matrix at input 

(EMF design) 


&. IMPLEMENTATION USING OUTPUT FEEDBACK 

Both the synthesis algorithms discussed herein lead to state feed- 
back control laws. Since in practice* 

not all the states of a system are measurable* an estimator design is 
required. Two issues that arise with estimation are the effect on tran- 
sient response and robustness. Two methods* which can be used to imple— 
ment the full-state control laws* that are worthy of consideration and 
will be pursued further in the future, are cited here. 

5 . 1. Robust Kalman Fi Iter 

This -procedure is as discussed in £7,83-and consists of parametri- 
cally increasing the process noise in the Kalman filter synthesis pro- 
cedure-till the full— state loop robustness is recovered. Properly 
implemented* this procedure does not increase the .dynamic order of the 
stick-response transfer functions. 

£• 2 Robust Output Observers 

This procedure is as discussed in £93 and consists of using 
observer theory to reconstruct system states from the available measur- 
ments. The advantages of this approach are that it leads to low order 
and guaranteed stable controllers and the sta te-f eedbac k system is fully 
reconstructed so there is no deterioration in performance and robust- 
ness. The procedure takes advantage of certain system structure* and 
therefore cannot always be applied* however. 
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6. Conclusions 

Some important conclusions based on the design examples addressed 
to date are as follows 

PEA 

The. direct eigenspace assignment technique has the advantage that 

t 

the resulting augmented system is of the same dynamic order as the 
open— loop system* which means that once the design is performed* the 
handling qualities evaluation of the augmented system can directly be 
carried out. Also* once the design requirements are properly mapped into 
the form of a desired eigenstructure*- the solution algorithm for obtain- 
ing the “closest possible" eigenstructure for the augmented system is 
quite straightforward. 

Some of the areas that require further investigation to make the 
DEA design technique more viable are 

a. Procedures for selecting the desired eigennstructure such that 
the choice reflects the control system design objectives. 

b. Procedures for improving the stability robustness of the feed- 
back design 

c. Prefilter design to properly blend the control inputs such 
that the desired dynamics are obtained from the pilot's stick 
input. 

EhF 

Though the EMF technique has the advantage that the flight 
control design requirements are easily mapped into the design pro— 
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cedure* its major disadvantages are that the resulting augmented 
system is of high order and the required bandwidth may be high. 
This means that system reduction is required in order to evaluate 
the handling qualities of the augmented system. This needs to be 
done for the control last synthesized herein. Finally* very high 
bandwidth control laws may be susceptible to model errors and 
require high deflection rates. 


/ 
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APPENDIX 

The macro for direct eigenspace assignment is listed in Table 
A. 1 and that for explicit model following in Table A. 2. Both the 
macros are written in the form of "user defined functions" and can 
be executed inside MATRIX^. Also, the macros are well documented 
with comment cards so that the user can understand the input-output 
requirements. 
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Table A. 1 Macro 'eigassgn' for DEA design 


// Cva. f b krl=eigassgn(a. b* nc. nsp, lame 1. vd* qd ) nDiriM-ai on 

// state feedback gains for eigenspace assignment GIIMAL PAGE !S 

// inputs are a* b* nc# nsp* lame 1» vd> qd 0F p OQR QUALITY 

// a — system matrix 

// b - system control distribution matrix 
// nc — no. of desired complex poles / 2 

// nsp - no. of complex poles for which the weighting on the eigenvector 
// erfor is different for real and imaginary parts 

’// lamcl — column vector of desired closed loop poles with nc+nr elements 
// where nr is the no. of. desired real poles 

// first nc elements are the desired complex poles with +ve imag. part 

// with the. nsp elements appearing first 

// vd — n by (nc+nr) matrix of desired eigenvectors with ith column corr. to 
// ith element of lamcl 

// qd — n by (nsp+nc+nr) matrix of weigthing on eigenvector errors* ith column 
// forms the diagonal elements of the weigthing matrix of ith eigenvector 

// the first 2#nsp columns correspond to the nsp eigenvectors with 

// weighting on the real and imaginary parts appearing consecutively 

// (the real part first) 

// outputs are* va* fbkr 

// va — matrix of acheivable eigenvectors 

// fbkr - feedback gain matrix (corr. to neg. feedback) 

// 

C temp 1 » temp23=s i 2 e ( lame 1 ) j nr=templ-nc» 
ns=2#nc+nr» 

// 

if nsp>0* Cns* nu3=s i 2 e ( b ) i bb=Cb O*ones (ns* nu ) * O#ones(ns.nu) bl*... 

for i = l : nsp* j=2*i-li qtmp = diag ( Cqd ( : * j >* qd ( : » j + 1 > 1 )i . . . 

vtmp = Cr eal ( vd ( : * i ) )* imag(vd(:* i))3* tmp l=r eal ( lame l(i))*eye(ns) -a* ... 

tmp2=imag ( lame 1 ( i ) )*ey e (ns )i mtmp = Ctmpl -tmp2» tmp2 tmplDi ... 

tmp =v tmp / *qtmp«-inv(mtmp )*bb*inv(bb '♦invCmtmp ) / *qtmp*inv(mtmp )*bb )i . . . 

vtmp = inv(mtmp )*bb*tmp '» tmp = tmp '* w ( : » j ) = tmp ( 1 : nu ) + jay *tmp (nu+1 : 2#nu) ; ... 

w(: * j + l) = conj(w(: » j))* va(: * j )=vtmp ( 1 : ns ) + jay*vtmp (ns + 1 : 2*ns ) ; ... 

va(: » j+1 )=conj ( va ( : . j ) )» . . . 

if i=nsp* clear qtmp bb tmp tmp 1 tmp2 mtmp vtmp nu* 

// 

tmp l=nc-nsp « 

if tmpl>0* for i=nsp + l:nc» j=2*i-l» 1 i = inv ( lame 1 ( i )*ey e (ns )-a ) *bi ... 

tmp=vd (: » i) '*diag(qd(: » nsp + i) )*li*inv(li '*diag(qd(: * nsp + i ) )*1 i ) i . . . 
w(:.j)=tmp'j w( : * j+1 )=con j (w( : '* j ) ) * va(: » j)=li*w(: * j)» va ( : * j+1 )=conj (va ( : * j ) ) 
if nr>0* for i = l:nrj j=2*nc + i» l=nc + i* 1 i = i nv ( lame 1 ( 1 )*ey e (ns ) -a )*bi ... 
tmp=vd (: » 1) / *diag(qd(: *nsp + l))*li*inv(li , *diag(qd(: *nsp + l))*li)j . . . 
w(:.j)=tmp'* va(:*j) = li*w(:,j)< 

// 

f b k=-w*inv ( va ) « f b kr=r ea 1 ( f b k ) * 
r etf 
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Table A. 2 Macro 'fcemf* for EMF design 


/ /Lew ah kr*saug,num] = fcemf(svism,tp,qez/r»hv) 

// flight control explicit model following algorithm 
// inputs are sv* sm> tp# qez* r< hv 
// sv — vehicle system; = Cav bv; cv O] 

// sm — model system; = Cam bm; cm 03 
// tp — pilot stick time constant 

// qez — weighting matrix for error .between vehicle and model outputs 
// r — control weighting matrix 

// hv — matrix of vehicle outputs for which transfer functions are desired 
// outputs are evalj kr* num 
// eval — closed loop eigenvalues 

// kr — state feedback gains corr. to neg. feedback 

// saug — augmented system; pilot stick as input and augmnetd vehicle outputs 
// num — numerator coeff. of augmented vehicle output transfer functionss 
// (ith row corr. to i th output) 

// 

Cnyj nvl=size(hv); 

C 1 1 , t2J=size(sv); nu=t2-nv; n z = t 1 -nv; C t3* nmi ]=size(sm); nm=nml-li 
a=Csv(l:nvi l:nv) 0*ones CnVi nml ) ; Ozones (nm, nv ) sm ( 1 : nm, 1 : nml )# 

0*ones < 1 1 nv+nm) , -1/tpD; b = Csv( 1 : nv, (nv+1 ) : t2); 0«ones(nml> null; 

c=Csv( (nv+1 ) : t 1 » 1 : nv ) » -sm(nml : t3» 1 : nm) 0#ones(nz ( 1)1; 
q=c '*qez*c; 

// 

Ceval, kr]=regulator(ai b. q, r ); 
temp=a-b*kr; eval = eig ( temp ); 

// 

t4=nm+nv; 

aaug = temp < 1: t4, 1 : t4); baug=temp < 1 : t4, t4+l ); 

caug = Chv 0#ones(ny. nm)]; saug=Caaug baug; caug 0*ones(nyi 1)1; 

Cnuni« den]=tform(saugi t4); 
retf 



